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Classification System for Sex Determination 

and Tissue Characterization 

BACKGROUND OF THE INVENTION 

TECHNICAL FIELD 

The invention relates to the determination of the sex of both human and 
animal subjects. More particularly, the invention relates to the instrumentation 
and method by which the sex and general tissue parameters of human and 
animal subjects can be determined noninvasively. 

DESCRIPTION OF THE PRIOR ART 

Near infrared (NIR) tissue spectroscopy is a promising noninvasive 
technology which bases measurements on the irradiation of a tissue site with 
NIR energy in the 700-2500 nanometer wavelength range. The energy is 
focused onto an area of the skin and propagates according to the scattering 
and absorption properties of the skin tissue. Therefore, the reflected or 
transmitted energy that escapes and is detected provides information about 
the tissue volume that is encountered. Specifically, the attenuation of the light 
energy at each wavelength is a function of the structural properties and 
chemical composition of the tissue. Tissue layers, each containing a unique 
heterogeneous particulate distribution, affect light absorbance through 
scattering. Chemical components such as water, protein, fat, and blood 
analytes absorb light proportionally to their concentration through unique 
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absorption profiles or signatures. The measurement of tissue properties, 
characteristics or composition is based on detecting the magnitude of light 
attenuation resulting from its respective scattering and/or absorption 
properties. 

5 

Blood Analyte Measurement 

While noninvasive measurement of blood analytes, such as blood glucose 
concentration, has been pursued through NIR spectroscopy, the reported 

1 0 success and product viability has been limited by the lack of a system for 
compensating for structural variations between individuals that produce 
dramatic changes in the optical properties of the tissue sample (for example 
see O. Khalil, Spectroscopic and clinical aspects of non-invasive glucose 
measurements, Clin. Chem., vol. 45, pp. 165-77 (1999) or J. Roe, B. Smoller. 

1 5 Bloodless Glucose Measurements, Critical Reviews in Therapeutic Drug 
Carrier Systems, vol. 15, no. 3, pp. 199-241 (1998). These differences are 
largely anatomical and provide distinct systematic spectral absorbance 
features or patterns that can be related directly to specific characteristics such 
as dermal thickness, protein levels, and percent body fat. While the 

20 absorbance features are repeatable by subject, over a population of subjects 
they produce confounding nonlinear spectral variation. Therefore, differences 
between subjects are a significant obstacle to the noninvasive measurement 
of blood analytes through NIR spectral absorbance. 

25 An apparatus and procedure for substantially reducing this problem by the 
classifying subjects according to major skin tissue characteristics prior to 
blood analyte prediction is described in S. Malin, T. Ruchti, An Intelligent 
System for Noninvasive Blood Analyte Prediction, U.S.- Patent Application 
Serial No. 09/359,191, filed July 22, 1999. The selected characteristics are 

30 representative of the actual tissue volume irradiated and the amount of the 
target analyte that is sampled. By grouping individuals according to the 
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similarity of spectral characteristics representing the tissue structure, the 
nonlinear variation described above is reduced and prediction of blood 
analytes becomes more accurate. 

5 In human subjects, significant differences related to sex have been 
discovered in the skin tissue. These differences include the thickness of the 
dermis (see C. Tan, B. Statham, R. Marks, P. Payne, Skin thickness 
measurement by pulsed ultrasound: its reproducibility, validation and 
variability British Journal of Dermatology, vol. 106, pp. 657-667, (1982) and 

1 0 J. Bliznak, T. Staple, Roentgenograph^ measurement of skin thickness in 
normal individuals, Radiology, vol. 116, pp. 55-60 (July 1975)), the amount of 
fat in subcutaneous tissue (see J. Durnin,.M. Rahaman, The assessment of 
the amount of fat in the human body from measurements of skinfold 
thickness, British Journal of Nutrition, vol. 21 (1967) and F. Johnston, 

15 Relationships between body composition and anthropometry Human 
Biology, Vol. 54, No. 2, pp. 221-245 (May 1982)) and skin collagen and 
density (see S. Shuster, M. Black, E. McVitie, The influence of age and sex 
on skin thickness, skin collagen and density, British Journal of Dermatology, 
vol. 93 (1975)). The determination of subject sex therefore provides an 

20 important indication of large systematic differences in the tissue structure and 
composition. 

Therefore, an automated method for the determination of the subject's sex 
provides valuable information relevant to subject classification and 
25 determination of key tissue properties for blood analyte measurement. 

Sex Determination of Animals 

The determination of the sex of animal species has commercial benefit in 
30 certain industries due to the replacement of a human expert by an accurate 
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and automated noninvasive device (see T. Miyakawa, O. Kato, Y. Koike, K. 
Matsunami , N. Sekiya, Fish sex discrimination equipment and method, 
United States Patent No. 5,013,906 (May 7, 1991); K. Suzuki, Apparatus for 
determining the sex of a chick, United States Patent No. 4,417,663 (Nov. 29, 
5. 1983); A. Frasch, R. Ugalde, Procedure for the sex determination of 
embryos in mammals especially applied to bovine embryos, United States 
Patent No. 5,578,449 (Nov. 26, 1996); and W. Cheng, C. Chen, C. Hu, C. 
Wang, K. Choo, Process for sexing cow embryos, United States Patent No. 
5,876,942 (March 2, 1999)). In Miyakawa et al supra, the sex of a fish is 

1 0 determined by examining the color of the genital gland area through visible 
light. In Suzuki et al supra, the sex of a chick is determined by examining the 
color of the anal region through the use of visible light. In Frasch et al and 
Cheng et al supra, methods for sexing cow embryos are detailed through a 
complex method of polymerase chain reactions. These methods are not 

1 5 extendable to human subjects or other mammals due to gross anatomical 
differences. Further, the methods are limited because they involve either the 
automated color detection of a particular often unexposed region of the animal 
or rely on measurements that are invasive or semi-invasive. Finally, none of 
the methods listed above use a near-infrared technology that penetrates the 

20 tissue to measure its internal properties. Therefore, a device for sex 
determination of animals needs to be developed that is accurate, noninvasive, 
automated, and general. 



Body Composition Determination 

25 

The automated and noninvasive determination of sex provides beneficial 
information related to the body composition of the subject. For example, in 
the determination of the lean-body mass of humans, the knowledge of the 
subject's sex is required prior to analysis of other anthropometric 
30 measurements (see V. Heyward, L. Stolarczyk,. Applied Body Composition 
Assessment, Human Kinetics, Champaign, IL (1996)). An automated and 
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noninvasive device for sex determination provides a critical component for a 
fully automated method of body composition analysis. 

It would be advantageous to provide a novel apparatus and related 
5 procedures for sex determination of human and animal subjects through NIR 
tissue spectroscopy that has particular benefit in several areas, including 
blood analyte prediction, animal sex determination, and body composition 
evaluation. 



10 

SUMMARY OF THE INVENTION 

The invention herein provides a novel method of sex determination for 
animals and humans based on near-infrared measurements of the skin tissue. 
15 In addition, the invention provides fundamental information regarding gross 
tissue characteristics and can be used for determination of systematic and 
relative differences in the thickness of the dermis and the amount of 
subcutaneous fat at the measurement site. 

20 The invention is a method for non-invasively determining the sex of human or 
animal subjects. The method uses a spectroscopic apparatus in conjunction 
with an optical interface to measure tissue properties and characteristics that 
are manifested spectrally and that vary systematically according to the 
subject's sex. 

25 

The procedure for sex determination involves a calibration model that is 
empirically derived from a set of exemplary samples consisting of NIR tissue 
measurements and the actual sex of a population of subjects. The model is a 
set of parameters and computer generated code that is implemented to 
30 predict the subject's sex. The prediction consists of a discrete sex 
determination (male or female) and one or more relative property magnitudes 
that reveal information regarding the tissue properties of the sampled tissue 
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volume. These properties include but are not limited to the thickness of the 
dermis, the collagen content, the skin density, and the amount of 
subcutaneous fat at the measurement site. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram of a sex classification instrument and sex 
classification procedure according to the invention; 

Figure 2 is a graph plotting a typical noninvasive absorbance spectrum; 

Figure 3 is a graph plotting NIR spectral measurements and demonstrating 
spectral shift between male and female subjects; 

Figure 4 is a block diagram showing a general procedure for determining 
individual sex based on the spectral shift of near infrared spectra according to 
the invention; 

Figure 5 provides histograms of the spectral minimum near 1270 nm 
separated by sex; 

Figure 6 provides histograms of the spectral minimum near 1650 nm 
separated by sex; 

Figure 7 provides histograms of the spectral minimum near 2200 nm 
separated by sex; 

Figure 8 is a block schematic diagram of a classification system for 
determining the sex of a subject based on an NIR measurement according to 
the invention; 



WO 01/52726 PCT/US00/34780 

Figure 9 provides plots of the scores from NIR spectral regions illustrating the 
separation of abstraction features related to sex; 

Figure 10 is a block schematic diagram of a procedure for sex determination 
5 through abstract features and linear classification according to the invention; 

Figure 11 is a graph showing loading of the third principal component 
compared to the absorbance spectrum of animal fat; 

1 0 Figure 12 is a graph showing NIR measurements from growth hormone mice 
using a spectrometer; 

Figure 13 is a plot showing separation by gender in principal component three 
scores; 

15 

Figure 14 is a plot showing a comparison of principal component three 
loading to a protein spectrum; 

Figure 15 is a plot showing a principal component three versus principal 
20 component one scores; 

Figure 16 is a block schematic diagram showing a procedure for sex 
determination on the basis of feature extraction through pure component 
absorbance spectra of water and fat according to the invention; 

25 

Figure 17 is a graph showing pure component absorbance spectra of water 
and triglycerides (fat); 

Figure 18 is a plot showing the extracted features related to the absorbance 
30 of fat and water separated by sex; and 
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Figure 19 is a block schematic diagram of a classification of subject sex 
based on NIR spectral measurements through optimal wavelength selection 
and a global model according to the invention. 

5 

DETAILED DESCRIPTION OF THE INVENTION 

The invention provides an apparatus for measuring the infrared absorption by 
tissue irradiated with near-infrared energy and a procedure of determining the 
1 0 subject's sex. 

Apparatus 

A block diagram of the integrated system is shown in Figure 1. The apparatus 
includes an energy source, a sensor element, an interface to the subject, a 

15 wavelength selection device and an analyzer. The source generates and 
transmits near-infrared energy in the wavelength range 700-2500 nanometers 
and consists of a device, such as an LED array or a quartz halogen lamp 12. 
The sensing elements include detector optics 15 and detectors 16 which are 
responsive to the targeted wavelengths. The method of wavelength 

20 separation includes a monochrometer 13, an interferometer or successive 
illumination through the elements of an LED array. The interface 14 to the 
subject 11 is a method of transmitting energy from the source to the target 
skin tissue measurement site and includes, for example a light pipe, fiber- 
optics, a lens system or a light directing mirror system. Energy is collected 

25 from the surrounding tissue areas in reflectance mode at an optimally 
determined distance(s) through the use of starring detectors or fiber optic 
probes. Alternately, energy is collected in a transmission mode through a skin 
flap, ear lobe, finger or other extremity. The collected light is converted to a 
voltage and sampled through an analog-to-digital converter 17 for analysis on 

30 a m icrop rocesso r based system . 

8 
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Gen ral Sex Classification Procedure 

The general procedure for sex determination based on the measured 
5 spectrum, shown in Figure 1, is implemented in a microprocessor that 
automatically receives the measurement information from the ADC 17. The 
principal components of the sex determination procedure include outlier 
detection 102, preprocessing 103, feature extraction 104, and classification 
105. The design of each procedure is performed on the basis of a calibration 
1 0 set of exemplary measurements. In this section, the general steps which are 
detailed in the subsequent Design and Implementation Sections are 
summarized. 



Measurement 

15 

The measurement is a spectrum denoted by the vector me^ of absorbance 
values pertaining to a set of N wavelengths AeSR" that span the near infrared 
(700 to 2500nm). A typical plot of m versus X is shown in Figure 2. 

2 0 Outlier Detection 

The outlier detection procedure is a method for detecting invalid 
measurements through spectral variations that result from problems in the 
instrument, poor sampling of the subject or a subject outside the calibration 
25 set. The preferred method for the detection of spectral outliers is through a 
principal components analysis and an analysis of the residuals (see H. 
Martens, T. Naes, Multivariate Calibration . John Wiley and Sons, New York, 
pp. 419 (1989)). First, the spectrum, m, is projected onto five eigenvectors, 
contained in the matrix o, that were previously developed through a principal 
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15 



components analysis (on a calibration set of exemplary absorbance spectra) 
and are stored in the computer system of the device. The calculation is given 
by 



xpc 0 = J^mo k (1 ) 



k=] 



and produces the 1 by 5 vector of scores, xpc 0 where o k is the kth column of 
the matrix o. The residual, q, is determined according to 



10 q =m-xpc 0 o T (2 ) 



and compared to three times the standard deviation of the expected residual 
(of the a calibration set). If greater, the sample is reported to be an outlier and 
the sex determination procedure is terminated. 



Preprocessing 



Preprocessing includes operations such as scaling, normalization, smoothing, 
derivatives, filtering, and other transformations that attenuate the noise and 
20 instrumental variation without affecting the signal of interest. The 
preprocessed measurement, xeSi N , is determined according to 



x = h(X,m) (3 ) 



25 where h:R Nx2 ->9? w is the preprocessing function. 
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Feature Extraction 

Feature extraction determines the salient characteristics of measurements 
that are relevant for sex determination. Feature extraction is any 
mathematical transformation that enhances a quality or aspect of the sample 
measurement for interpretation. The purpose of feature extraction is to 
represent concisely and enhance the properties and characteristics of the 
tissue measurement site for sex determination. In addition, the features 
provide significant information of the tissue properties they represent and can 
be used for alternate purposes, such as diagnostics or system optimization. 

The features are represented in a vector, ze9t M that is determined from the 
preprocessed measurement through 



* = (4 ) 

where t 9i w ->9l M is a mapping from the measurement space to the feature 
space. Decomposing f(»)yields specific transformations, f{*)\ 9t w -»9t M , for 
determining a specific feature. The dimension, M h indicates whether the ith 
feature is a scalar or a vector and the aggregation of all features is the vector 
z. When a feature is represented as a vector or a pattern, it exhibits a certain 
structure indicative of an underlying physical phenomenon. 

The individual features are divided into two categories: abstract and simple. 
Abstract features do not necessarily have a specific interpretation related to 
the physical system. Specifically, the scores of a principal component 
analysis are useful features although their physical interpretation is not always 
known (see H. Martens, T. Naes, supra.). For example, the utility of the 
principal component analysis is related to the nature of the tissue absorbance 
spectrum. The most significant variation is generally related to the structure 
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which varies systematically with sex. Therefore, the scores of the principal 
components analysis provides useful information for sex determination and 
constitute a valuable set of features. 

5 Simple features are derived from an a priori understanding of the sample and 
can be related directly to a physical phenomenon. For example, the thickness 
of the dermis or subcutaneous layer, described previously as varying 
systematically with sex, and result in specific spectral manifestations. These 
spectral variations are extracted and enhanced and serve as both a feature 
1 0 for sex determination and a measurement of their respective tissue properties. 

In the general case, the full spectrum can be passed to the classification 
system for sex determination. However, the presently preferred embodiment 
of the invention provides three specific methods of feature extraction that 
15 exhibit superior classification performance and measurements of other 
relevant tissue properties: 

1 . The scores from factor analysis; 

20 2. Location of critical points on the measured spectrum; and 

3. Relative absorption of water and fat. 

25 The detailed implementation of the procedure for extracting these features on 
the basis of a calibration set is provided in the next section. 

Classification 

30 The determination of the subject's sex on the basis of the extracted features 
is performed through a classification step which involves a mapping and a 
decision. The mapping step is given by 

1 3 
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L = /go 



(5 ) 



where L is a scalar that can be used to measure the distance for the 
categories of male and female. Generally, two values, L male and L femgle , 
associated with the representative or mean value of L for the male and female 
categories respectively are predefined and the class assignment is based on 
the closeness of L to L maie and L lemale . For example, the distance of L to 
previously established class means classes can be measured by 



The mapping and decision limits are determined from a calibration set of 
exemplary features and corresponding sexes through a classification 
calibration procedure. Existing methods include linear discriminant analysis 
(see R. Duda, P. Hart, Pattern Classification and Scene Analysis. John Wiley 
and Sons, New York (1973)), SIMCA (see S. Wold, M. Sjostrom, SIMCA: A 
method for analyzing chemical data in terms of similarity and analogy, 
Chemometrics: Theory and Application, ed. B.R. Kowalski, ACS Symposium 
Series, vol. 52 (1977)), k nearest-neighbor (see R. Duda, P. Hart, supra.), 
fuzzy classification (see J. Bezdek, S. Pal, eds., Fuzzv Models for Pattern 
Recognition, IEEE Press, Piscataway, NJ (1992) and J. Keller, M. Gray, J. 




(6 ) 



the decision is made as follows 



if d ma , g <d fema i e then male 
i f oUfc^cffemate then female 
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Givens, A Fuzzy K nearest Neighbor Algorithm, IEEE Transactions on 
Systems, Man, and Cybernetics, Vol. SMC-15, No. 4, pp. 580-585 
(July/August, 1985)), and various forms of artificial neural networks (see Y. 
Pao, Adaptive Pattern Recognition and Neural Networks . Addison-Wesley 
5 Publishing Company, Inc., Reading, MA (1989)). 

Implementation Details 

This section discloses four specific procedures for sex determination. The 
1 0 structure of the procedures are based on a priori knowledge of the systematic 
variation of the skin structure, i.e. the variation of skin thickness, amount of 
subcutaneous fat and collagen structure. However, the parameters of each 
procedure, such as the eigenvectors for outlier detection, are determined on 
the basis of a experimental data set providing exemplary information. 

15 

Experimental Data Set 

A study was performed to generate calibration and validation data for the four 
procedures disclosed subsequently. Human subjects (266) of diverse age, 

20 sex, and ethnicity were recruited at a local health care facility and detailed 
demographic information about each participant was recorded. Four replicate 
absorbance spectra were measured on each subject's forearm with the 
spectrometer described as the preferred embodiment The number of 
samples per participant was limited to one. The subjects were separated at 

25 random into calibration and test sets. The calibration set is employed to 
construct the models, mappings, and parameters for each of the procedures 
described below. The validation set was used to test the performance of each 
procedure for sex determination. The total set of spectra and demographic 
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information shall be referred to as the "Experimental Data Set" in the 
remainder of the text. 

While this is a specific experiment aimed at the determination of a suitable set 
for calibrating the sex determination apparatus, one can readily appreciate 
that for different subjects and for different target performance levels other 
experiments with more or less subjects would be performed. 



Method 1 - Sex Determination through Spectral Shift 

The first method predicts sex based on the spectral shift observed in the 
measured NIR spectrum. The basis for the method is that systematic sex 
related differences in the characteristics and properties of the subcutaneous 
fat and dermis layers causes systematic variation of the absorbance bands of 
specific analytes, such as fat. However, other background analytes tend to 
remain constant. The result of the variation in one absorbance band among a 
particular background is manifested through an apparent shift in the peak and 
valley locations. 

For example, the spectra of the Experimental Data Set was separated 
according to sex and averaged. The averaged spectra are depicted in Figure 
3 and show a pronounced difference between the sexes. Specifically, 
particular peaks and valleys appear shifted between the two sexes as marked 
on the figure. The wavelength position of the peaks and valleys (the critical 
points) provide a set of features suitable for determination of sex. 

The general procedure for utilizing these features is depicted in Figure 4. 
First, NIR absorbance spectra are measured 40, such as those shown in 
Figures 2 and 3. The measured spectrum is preprocessed 41 via a 15-point 
Savisky-Golay smoother in the form of a finite impulse response filter (see A. 
Savitzky, M. Golay, Smoothing and Differentiation of Data by Simplified Least 
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Squares Procedures, Anal. Chem., vol. 36, no. 8, pp. 1627-1639 (1964)). 
The apparent shift of the measurement at the critical points is measured 42 
with reference to a processed standard absorbance spectrum 43. The 
standard absorbance spectrum can be any single absorbance measurement 
or the average of a set of measurements. In the preferred embodiment of the 
invention, the standard absorbance spectrum is the mean of all spectra in the 
calibration set. 

A statistical model 45 is used to identify critical points that are applied using a 
decision rule 44. The location of the critical points is determined over a finite 
window in the vicinity of the known critical points (for example, a 10 nm 
window is used). Locations for critical points that vary according to subject 
sex include: 1145, 1195, 1230, 1270, 1650, 2160, 2200, 2250, 2305 and 
2350 nm. However, the method is general and can be applied to other 
spectral regions. In addition, other critical points can be easily selected from 
the first or second derivative of the absorbance spectrum. 

The wavelength position of each critical point is determined by fitting a 
second-order polynomial function to the measured spectrum in the vicinity of 
the known critical point locations, determining the derivative of the polynomial, 
and calculating the root. For example, given the polynomial 



y = a + bx + cx 2 ( 7 ) 

where x is the wavelength and y is the estimated absorbance in the vicinity of 
a critical point. The parameters a, b and c are calculated through standard 
regression techniques between xand y, the measured absorbance spectrum. 
The estimated peak location, p, is given by 
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(8 ) 



Alternate methods of determining the critical point can also be used, such as 
zero-crossing algorithms, gradient search algorithms, or methods based on a 
5 cross-correlation function between the standard and measurement spectra. 

Given a calibration set of spectra, the spectral features (critical points) are 
extracted as described above and linear discriminant analysis or Mahalanobis 
distance is used to determine the classification model 44 (see R. Duda, P. 
10 Hart, supra.). For example, the Mahalanobis distance of a set of critical 
points associated with a measured spectrum is given by 

d moh =(P- Pmah) S ' X (P ~ Pmale) ^ g 

d f m ak = (P " P female)^ ' (P ~~ P female) 

15 were p is the vector of critical points, p mate is the mean of the critical point 
locations over all males in the calibration set, p female is the mean of the critical 
point locations over all females in the calibration set, S is the covariance of 
the critical points over the calibration set and d ma)e and d femal0 correspond to 
the closeness of the critical point to the male and female class represented by 

20 the calibration set. The same critical points are calculated for future spectral 
measurements and Equation 9 is applied to determine the sex of the subject. 

For classification using a single critical point, the squared statistical distance 
is calculated and used to determine the sex of the subject through the 
25 calculation 
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j2 
male 



d 2 

female 



f V 

P ~ Pn,ale 
V S male J 

P P female 
^ S female J 



(10 ) 



where p is the location of the critical point, p male and s male are the mean and 
standard deviation of the critical point locations over all males in the 
5 calibration set, p femalQ and s female are the mean and standard deviation of the 
critical point locations over all females in the calibration set and d ma)e and d female 
correspond to the closeness of the critical point to the male and female class 
represented by the calibration set. 

10 Given the distances calculated in Equations (9) or (10), the sex is determined 
46 (Figure 4) through 



15 



if ^na/^^fema/e then male 
if d ma ,e>d fema}e then female 

The same critical points are calculated for future spectral measurements and 
Equation 9 or 10 is applied with these decision rules to determine the sex of 
the subject. 



20 As an example, the procedure for sex determination was applied to the 
calibration and test sets comprising the Experimental Data Set. First, the 
critical points were determined at the following wavelengths 1270, 1650, and 
2200 nm. The calculated critical points for the each sample of the 
Experimental Data Set are plotted via histograms in Figures 5-7. The figures 

25 show that the locations of the critical points are normally distributed about a 
mean level. However, this mean is systematically different between males 
and females. The determination of sex is possible through the process 
outlined above and produced validation set results on the Experimental Data 
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Set with a sex determination accuracy of 85%. The use of additional critical 
points further improves the performance, while fewer critical points afford 
simpler procedures but degrade performance. 

5 Furthermore, the various critical points shown provide an indication of the 
thickness of the dermis, subcutaneous tissue and the properties of the dermis. 
For example, the points at 1270 nm vary primarily according to the thickness 
of the subcutaneous fat due to the deeper penetration of the NIR energy. The 
critical point at 1 650 can be used to provide relative information about the 
1 0 thickness of the dermis and the critical point at 2200 nm is likely related to the 
properties of collagen fibers in the dermis. This information can be used to 
further classify the subjects for blood analyte prediction, analysis of body 
composition or determination of environmental influences on skin properties. 

1 5 Method 2 - Sex Determination through Abstract Feature Extraction 

The second method, illustrated in Figure 8, is distinguished by the use of 
factor analysis to develop a set of abstract features capable of representing 
the spectral variation related to sex. The measurement is a NIR absorbance 

20 spectrum 80 similar to that shown on Figure 2. The spectrum is sub-divided 
into one or more regions according to wavelength (wavelength selection) and 
is preprocessed and normalized 81 to enhance spectral variation related to 
sex. The measurements are projected 82 onto one or more sets of previously 
determined factors (eigenvectors) 83 to determine the scores. The scores are 

25 the extracted features and are subjected to a classification procedure 84-87, 
such as linear discriminant analysis, SIMCA, k nearest-neighbor and various 
forms of artificial neural networks to predict the sex of the subject. 

Example 1 - Human Subject Sex Determination 

30 As an example the Experimental Data Set was analyzed to demonstrate the 
utility of abstract features. First, wavelength selection was applied to reduce 

20 
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the spectra to the 1 100-1400 nm range. An 1 1 -point Savisky-Golay smoother 
was applied and the set of spectra were subjected to multiplicative scatter 
correction, as described below. A principal components analysis was 
performed and the scores of the first three eigenvectors were calculated and 
5 are plotted in Figure 9 (90, 92, 94) according to sex. The scores, 
representing variation in the spectra, show a pronounced systematic 
separation according to subject's sex. 

The preferred embodiment of the sex determination (through the abstract 
1 0 feature extraction) procedure is shown in Figure 1 0 and involves spectral 
preprocessing 110-1 12, decomposition through principal components analysis 
113, and classification through linear discriminant analysis 1 14-117. 

First, the absorbance spectrum, m, is provided from the outlier detection 
15 system 110. Wavelength selection 112 is applied to truncate the spectral 
range to regions with significant absorption due to fat in adipose tissue (1100 
to 1400 nm). The spectrum is also processed through multivariate scatter 
correction 1 1 1 (see P. Geladi, D. McDougal, H. Martens, Linearization and 
Scatter-Correction for Near-Infrared Reflectance Spectra of Meat, Applied 
20 Spectroscopy, vol. 39, pp. 491-500 (1985)) through a rotation that fits it to the 
expected or reference spectrum, m , determined from the calibration set. 
First, the spectrum is fit via linear regression according to 

m = a+bm + e / -\-\ 

25 

where a and b are the slope and intercept and e is the error in the fit. The 
spectrum is then corrected through 
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where x is the processed absorbance spectrum. The processed spectrum is 
projected onto the eigenvectors, p h that were previously developed through a 
principal components analysis (on the calibration set) 113. The calculation, 
5 shown in Figure 1 0, produces the 1 by N vector of scores, xpc. 

A discriminant function is applied to classify the subjects on the basis of the 
first M scores {M=5 is this application). The scores are rotated through a 
cross product with the discriminant , w, as depicted in Figure 10 to produce 
10 the scalar, L 1 14. This result is compared to I, the center between the two 
classes 115. If L > L then the subject is classified as a female 117. If not, the 
spectrum is classified as beloning to a male 116. 

This system was implemented with five factors selected and the mean 
1 5 spectrum for MSC, the eigenvectors (p) and discriminant function (w) shown, 
were calculated using the calibration set. Application of the procedure to 
determine the sex of the validation set resulted in a 90% correct sex 
determination accuracy. 

20 The method of classification described above can be replaced by other more 
suitable techniques depending on the resource and performance 
requirements. For example, the Fuzzy-K nearest neighbor classification 
algorithm reported by Keller (see J. Keller, M. Gray, J. Givens, supra.) was 
applied with 15 neighbors in place of the discriminant analysis discussed 

25 above. The method resulted in a validation set accuracy of 93%. 

Further, the scores of the abstract factors provide relative information 

i 

regarding the thickness of the subcutaneous fat, the dermis, and the level of 
collagen density in the dermis. For example, the factor corresponding to the 
30 third principal component of Figure 9 is plotted in Figure 11 along with the 
pure component spectrum of animal (bovine) fat. The correspondence of the 
major fat absorbance bands between the two features indicates that the 
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scores of the third factor represent the absorbance due to fat manifested in 
the measured spectra. The differences between the factor and the fat 
spectrum are caused by the covariation of other interferents and differences in 
the mode of measurement. 

5 Example 2 - Sex Determination of Mice 

Sex Determination through Abstract Feature Extraction was implemented for 
the determination of the sex of genetically altered growth hormone mouse. 
Growth hormone transgenic mice have foreign growth hormone genes 

10 introduced into their genome under the control of a promoter. The 
introduction of foreign growth hormone genes into their genome causes the 
animal continually to produce excessive amounts of growth hormone (see R. 
Wanke, E. Wolf, W. Hermanns, S. Folger, T. Buchmuller, G. Brem, The GH- 
Transgenic Mouse as an Experimental Model for Growth Research: Clinical 

1 5 and Pathological Studies, Hormone Research, vol. 37, pp. 74-87(1992)). 

Spectral data was collected on a diffuse reflection spectrometer from fourteen 
growth hormone mice in the wavelength range 1100-2400 nm. The belly of 
the mouse was scanned and the intensity spectra were converted to 
20 absorbance as described previously. A plot of the resulting absorbances by 
sex revealed an overall higher absorbance across all wavelengths on female 
mice as compared to male mice (Figure 12). Based on the plot, wavelength 
selection was performed and the spectra were narrowed to 1500 to 2450nm. 

25 Principal component analysis (PCA) was performed on the selected 
wavelength region (1500 to 2450nm) of the growth hormone mouse data. 
The first three principal components were found to be significant through 
inspection of the factor scores. Principal component three scores separate 
into two groups according to sex (Figure 13). Principal component three 

30 loadings resemble protein spectrum (Figure 14). This indicates the difference 

23 



WO 01/52726 



PCT/US00/34780 



between the mouse sexes are related to variations in the protein content of 
the mouse tissue. The mice cluster into two groups corresponding to sex by 
plotting principal component one scores against principal component three 
scores (Figure 15). 



10 



Linear discriminant analysis, described above, was applied to the first three 
principal components scores. Sex determination performance was evaluated 
through cross-validation using the "leave-one-ouf strategy and calculating the 
percent prediction error. The cross-validation procedure is used iteratively to 
predict the mouse sex by using all other samples to construct the calibration 
model (see H. Martens, T. Naes, supra.). Using this technique, the growth 
hormone mice were correctly classified with 88% test accuracy (see Table t). 



15 



Table 1 : Gender Classification Results 



Classification Variable 




GH Mice 


Number of Classes 




2 


Number of Correct Classifications 




32 


Number of Misclassifications 




4 


Percent Correctly Classified 




89 


Percent Error 




11 


Results from 


Classification 




Class Misclassified 


Correct 


Total 


1 2 


14 


16 


2 2 


18 


20 



20 



25 



30 



Therefore, the growth hormone mice can be classified according to gender by 
differences in the structural and chemical composition of the protein in their 
tissue. Further, the systematic differences in sex are modeled by absorption 
features corresponding to protein. Therefore, the abstract features provide 
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quantitative and relative information about the significant properties of mice 
tissue related to sex. 



Method 3 - Sex Determination through Projection of a Basis Set 

5 

The third method involves the projection of a basis set of pure component 
spectra onto the measured noninvasive absorbance spectrum to generate key 
features. The procedure, depicted in Figure 16, involves the preprocessing 
and wavelength selection of the measured absorbance spectrum 160. A 

10 basis set of pure component spectra 161 is projected onto the measured 
spectrum 162, as shown in the figure. The basis set is comprised of the key 
constituents, such as water, fat, and protein, which are related to the 
systematic spectral variation caused by the sex of the subject. In the figure, 
the absorbance spectra of water and fat are used to create the 2 by M vector, 

15 p, comprising the basis set. However, additional analyte spectra can be 
added for additional information or accuracy, for example protein. 

The basis set is fit to the measured spectrum, as shown in the figure, and 
produces the 1 by 2 vector of weights, m. The weights are representative of 
20 the absorbance due to water or fat and provide relative information about the 
thickness of the dermis (in the case of the magnitude of water) and the 
amount of subcutaneous fat (in the case of the magnitude of fat) 163. This 
information can be used for characterizing the subject, body composition 
assessment, or for diagnostic purposes 164. 

25 

A discriminant function 1 65 is applied to classify the subjects based on the 
two features contained in m through the equation shown in the figure to 
produce the scalar, L This result is compared to l, the center between the 
two classes 166. If l>L then the subject is classified as a female 167. If 
3 0 not, the spectrum is classified as beloning to a male 1 68. 

25 
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As an example, the pure component spectra shown in Figure 1 7 of fat and 
water were normalized to unit length and applied to the Experimental Data Set 
according to the procedure in Figure 16. The wavelength regions were 
limited to 1 100-1350 nm and the relative magnitude of water and fat that was 
5 produced by projecting the basis set onto each absorbance spectrum in 
shown in Figure 18. As Figure 18 indicates, the absorbance of water is higher 
in males than in females. Since the absorbance in the dermis is primarily due 
to the concentration of water, and the concentration of water is within a limited 
range between subjects, the relative magnitude of water indicates the 

1 0 pathlength through the dermis. This is consistent with published reports that 
the dermis of males is thicker than the dermis of females (see C. Tan, B. 
Statham, R. Marks, P. Payne, supra.). In addition, the magnitude of fat is 
systematically higher in females than males indicating greater absorbance 
due to the fat in subcutaneous tissue. Further, this general procedure can be 

1 5 used to extract any specific characteristic that distinguishable by a unique or 
dominant absorber. 

Using the calibration set, the linear mapping of Figure 16 was developed via 
linear discriminant analysis (see R. Duda, P. Hart, supra.) and produced the 
20 following weight vector 



w= [-0.599 0.9982] . ( ■) 

25 From the calibration set, the mean value for L was found to be 1 .51. Using 
these parameters in conjunction with Figure 18 produced a sex prediction 
accuracy in the validation set of 87%. 

Method 4 
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The final method disclosed (see Fig. 19) uses a general linear or nonlinear 
mapping to directly determine the sex of the subject. The NIR spectral 
measurement 190 is subjected to preprocessing and normalization 191 to 
reduce background interference and enhance variation related to sex. The 
key wavelength regions related to sex are extracted 192, which include the 
first and second overtone regions (1100-1800 nm), and provided to a 
classification model 193. The model is a linear or nonlinear mapping which 
produces the subject sex classification 194 directly. 



0 Although the invention is described herein with reference to the preferred 
embodiment, one skilled in the art will readily appreciate that other 
applications may be substituted for those set forth herein without departing 
from the spirit and scope of the present invention. Accordingly, the invention 
should only be limited by the Claims included below. 
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CLAIMS 

1. A method for determining the sex of animal and/or human subjects 
5 based on measurements of a skin tissue, comprising the steps of: 

measuring tissue properties and characteristics that are manifested 
spectrally and that vary systematically according to said subject's sex; 

applying a calibration model that is empirically derived from a set of 
exemplary samples consisting of tissue measurements and the actual sex of a 
1 0 population of subjects; and 

determining said subject's sex. 

2. The method of Claim 1, wherein said measuring step comprises the 
step of: 

15 using a spectroscopic apparatus in conjunction with an optical 

interface. 

3. The method of Claim 1, wherein said measurements are made in a 
spectrum denoted by the vector me9i w of absorbance values pertaining to a 

20 set of N wavelengths XeW" that span the near infrared (700 to 2500nm). 

4. The method of Claim 1, wherein said method further comprises the 
step of: 

determining tissue properties of a sampled tissue volume. 

25 

5. The method of Claim 4, wherein said properties comprise any of 
thickness of the dermis, collagen content, skin density, and the amount of 
subcutaneous fat at a measurement site. 

30 6. The method of Claim 1, wherein said determining step further 
comprises any of the steps of: 
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performing an outlier detection step to detect invalid measurements 
through spectral variations that result from any of problems in the instrument, 
poor sampling of the subject, and a subject outside the calibration set; 

performing a preprocessing step; 
5 performing a feature extraction step to represent concisely and 

enhance the properties and characteristics of a tissue measurement site for 
sex determination; and 

performing a classification step. 

1 0 7. The method of Claim 6, wherein said outlier detection step comprises 
the step of: 

performing a principal components analysis and an analysis of the 
residuals. 

15 8. The method of Claim 6, wherein said preprocessing step comprises 
any of the steps of: 

performing scaling, normalization, smoothing, derivatives, filtering, and 
other transformations that attenuate the noise and instrumental variation 
without affecting the signal of interest. 

20 

9. The method of Claim 6, wherein said feature extraction step comprises 
the step of: 

applying a mathematical transformation to enhance a quality and/or 
aspect of a sample measurement for interpretation. 

25 

10. The method of Claim 9, wherein said mathematical transformation 
comprises any of: 

a scores from factor analysis; 

location of critical points on a measured spectrum; and 
30 a determination of relative absorption of water and fat. 
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11. The method of Claim 6, wherein said classification step comprises the 
steps of: 

performing a mapping; and 
making a decision. 

12. The method of Claim 11, wherein mapping and decision limits are 
determined from a calibration set of exemplary features and corresponding 
sexes through a classification calibration procedure. 



10 13. The method of Claim 12, wherein said calibration procedure comprises 
any of linear discriminant analysis, k nearest-neighbor, fuzzy classification, 
and the use of artificial neural networks. 



14. A method for determining the sex of animal and/or human subjects 
1 5 based on measurements of a skin tissue, comprising the steps of: 

measuring an NIR spectrum of said subject's skin tissue; and 
predicting sex based on spectral shift observed in said measured NIR 
spectrum; 

wherein systematic sex related differences in characteristics and 
20 properties of said subject's subcutaneous fat and dermis layers causes 
systematic variation of absorbance bands of specific analytes while other 
background analytes tend to remain constant; and 

wherein variation in one absorbance band among a particular 
background is manifested through an apparent shift in peak and valley 
25 locations. 



1 5. The method of Claim 1 4, further comprising the step of: 
preprocessing said measured spectrum via a 15-point Savisky-Golay 

smoother in the form of a finite impulse response filter. 

1 6. The method of Claim 1 5, further comprising the step of: 
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measuring an apparent shift at a critical points with reference to a 
processed standard absorbance spectrum. 

17. The method of Claim 16, further comprising the step of: 
5 using a statistical model to identify critical points that are applied using 

a decision rule, wherein a location of said critical points is determined over a 
finite window in the vicinity of known critical points and wherein locations for 
critical points that vary according to subject sex. 

10 18. The method of Claim 1 7, further comprising the steps of: 

determining a wavelength position of each critical point by fitting a 
second-order polynomial function to the measured spectrum in the vicinity of 
said known critical point locations; 

determining a derivative of said polynomial; and 
1 5 calculating a root. 

1 9. The method of Claim 1 8, further comprising the step of: 
extracting spectral features. 

20 20. The method of Claim 1 9, further comprising the step of: 

using linear discriminant analysis or Mahalanobis distance to determine 
a classification model. 

21 . A method for determining the sex of animal and/or human subjects 
25 based on measurements of a skin tissue, comprising the steps of: 

measuring an NIR spectrum of said subject's skin tissue; and 
using factor analysis to develop a set of abstract features capable of 

representing the spectral variation related to sex. 

3 0 22. The method of Claim 21 , further comprising the steps of: 

sub-dividing said spectrum into one or more regions according to 
wavelength (wavelength selection); 
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preprocessing said regions; and 

normalizing said regions to enhance spectral variation related to sex. 

23. The method of Claim 22, further comprising the step of: 

projecting measurements of said regions onto one or more sets of 
previously determined factors to determine scores, wherein said scores are 
extracted features. 

24. The method of Claim 23, further comprising the step of: 

subjecting said scores to a classification procedure that, optionally, 
comprises any of linear discriminant analysis, SIMCA, k nearest-neighbor, 
and a form of artificial neural networks, to predict the sex of the subject. 

25. The method of Claim 14, further comprising the steps of: 
spectral preprocessing; 

decomposition through principal components analysis; and 
classification through linear discriminant analysis. 

26. A method for determining the sex of animal and/or human subjects 
based on measurements of a skin tissue, comprising the steps of: 

measuring an NIR spectrum of said subject's skin tissue; and 
projecting a basis set of pure component spectra onto a measured 
noninvasive absorbance spectrum to generate key features. 

27. The method of Claim 26, further comprising the step of: 
preprocessing and wavelength selection of said measured spectrum. 

28. The method of Claim 26, wherein said basis set is comprised of key 
constituents which are related to the systematic spectral variation caused by 
the sex of a subject. 
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29. The method of Claim 26, further comprising the step of: 
applying a discriminant function to classify subjects. 

30. A method for determining the sex of animal and/or human subjects 
5 based on measurements of a skin tissue, comprising the steps of: 

measuring an NIR spectrum of said subject's skin tissue; and 
using either of a general linear or nonlinear mapping to directly 
determine the sex of a subject. 

0 31 . The method of Claim 30, further comprising the step of: 

subjecting said NIR spectral measurement to preprocessing and 
normalization to reduce background interference and enhance variation 
related to sex. 

5 32. An apparatus for determining the sex of animal and/or human subjects 
based on measurements of skin tissue, comprising: 

means for measuring tissue properties and characteristics that are 
manifested spectrally and that vary systematically according to said subject's 
sex; 

means for applying a calibration model that is empirically derived from 
a set of exemplary samples consisting of tissue measurements and the actual 
sex of a population of subjects; and 

means for determining said subject's sex. 

33. An apparatus for determining the tissue properties of animal and/or 
human subjects based on measurements of skin tissue, comprising: 

means for measuring tissue properties and characteristics that are 
manifested spectrally and that vary systematically according to a selected 
characteristic of said subject; 

means for applying a calibration model that is empirically derived from 
a set of exemplary samples consisting of tissue measurements and the actual 
characteristic of a population of subjects; and 
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means for determining said characteristic of said subject. 



WO 01/52726 



PCT/US00/34780 



11/19 



Loadings Compared to Normalized Animal Fat 
0.4t 




-1-1 , , 1 

1000 1500 2000 2500 

Wavelength (nm) 

FIG. 11 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCTYUS00/34780 



1/19 



10 



11 



1 



16 



Instrument 
Components 




Monochrometer 



I 







Sub ject 


Sub ject 




Interface 




Module 


14- 
15- 


/ \ 




Detector Optics 



r 



Detectors and 
Analog 

Electronics 





S6x 

Classification 
Procedure 



100 



Absorbance 
Calculation 



T 



Outlier 
Detection 



T 



Preprocessing 



i 



Feature 
Extraction 



101 



102 



103 



104 




FIG. 1 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



2/19 




Wavelength(nm) 

FIG. 2 

(PRIOR ART) 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



3/19 




1000 1200 1400 1600 1800 2000 2200 2400 



FIG. 3 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/USOO/34780 



4/19 



40 



\ 



NIR Spectral 
Measurement 



I 



m 



43-^ 


4K 


Preprocessing 
and 

Normalization 


x 


Standard 
Noninvasive 
Spectrum 




Detection of 




Spectral Shift 






Statistical 




Decision Rule 


Model 




45-1 


U 




Sub feet Sex 
Classification 



42 



44 



46 



FIG. 4 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



5/19 



CO 

c 

D 
O 

CO 



CD 



Women 




1264 1266 1268 1270 1272 1274 1276 1278 




t test results: 
Different Means 
99.9% Confidence Int 
Prob. of Mistake=0 



mm 



1264 1266 



1268 1270 1272 1274 1276 
Location of Minimum (nm) 



1278 



FIG. 5 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



6/19 



PCT/US00/34780 



Women 




1642 1644 1646 1648 1650 1652 1654 1656 1658 1660 



Location of Minimum (nm) 

FIG. 6 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



7/19 




2192 2194 2196 2198 2200 2202 2204 2206 



Location of Minimum (nm) 

FIG. 7 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



8/19 



83 



1 



Factors 



86 



Statistical 
Model 



NIR Spectral 
Measurement 



| m 



Preprocessing 
and 

Normalization 



. Factor Analysis 



Sub ject 
Classification 



80 



■81 



82 




85 



Subject Sex 
Classification 



87 



FIG. 8 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



9/19 




—i LXt i i 1 i i i i I CM 




Scores of PC2 Scores of PC3 



SUBSTITUTE SHEET (RULE 26) 



PCT/US00/34780 



10/19 



NIR Spectrum 



110 



I 



m 



Multiplicative 

Scatter 
Correction 

i 

Wavelength 
Selection 



111 



112 



PCA 
N 

xpcj= I P k jx k 



k=1 



I 



Linear Mapping 
M 

L = I w k xpc k 



113 




Female 



114 



Male 



116 



FIG. 10 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/USOO/34780 



12/19 



CD 
O 

c 

€1-4 

o 

(/) 

< 

1.2 



Male and Female Growth Hormone Mice Spectra 



Males 



Females 




1500 1600 1700 



1800 1900 2000 2100 2200 2300 
Wavelengths 



2400 



FIG. 12 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



13/19 



Principal Component 3: Male and Female Growth Hormone Mice 
0.8 



0) 0.6 
O 

o _ 

^ 0.4 
0.2 

0} 



. O 



L°0°°° 0 ° * 
O ^ 



o 
o 



0 



-0.2 



D 

a. 

o -0.4 

^ -0.6 
-0.8 



O 



** * 



O Males 
* Females 



10 



15 20 25 30 
Samples 



35 40 



FIG. 13 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



14/19 



Principal Component 3 Loadings Compared to a Pure Protein Spectrum 

A 



1.5 



PC 3 Loadings 
Protein 



/ a\ 




1500 1600 1700 1800 1900 2000 2100 2200 2300 2400 

Wave lengths 



FIG. 14 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



15/19 



Growth Hormone Mices: PC 3 Scores vs. PC 1 Scores 



0.8 r 



CO 
Q 

o 
o 

(/) 



c 



0.6 



0.4 



0.2- 



0 



o 

Q. 

o -0.2 
o 



a -o.4 

o 

i= -0.6 



-0.8 



30 



o o 



Oq 0 
* O 



O Males 
♦ Females 



* * 



35 40 45 50 

Principal Component 1 Scores 



55 



FIG. 15 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



16/19 



PCT/US00/34780 



164 



\ 



L 



161 



Relative information about the 
thickness of the dermis and 
amount of subcutaneous fat. 



160 



NIR Spectral 
Measurements 



Pure 




Component 


P 


Water and Fat 




Spectra 





X 



Pro |ection 
M>P T ]~ 1 P X 



162 



m 





Relative 

Magnitude of 

Water and Fat 
Absorbance 










Linear Mapping 

N 

L= S w k m k 

k=1 




IL s-1 



■163 



165 




Female 



FIG. 16 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCTAJS00/34780 



17/19 



Pure Component Absorbance Spectra of Triglycerides and Water 



0.1 
0.09 

0.08 

0.07 



g 0.06 

C 

-e cos 

o 

a o.o4h 



o 



Normalized Absorbance of Water 
Normalized Absorbance of Triglycerides 




1200 1400 1600 1800 2000 2200 2400 

Wavelength 



FIG. 17 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/US00/34780 



18/19 




1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 

Relative Absorbance of Water 



FIG. 18 



SUBSTITUTE SHEET (RULE 26) 



WO 01/52726 



PCT/USOO/34780 



19/19 



N1R Spectral 
Measurement 



] m 



Preprocessing 
and 

Normalization 



Optimal 
Wavelength 
Selection 



Subject Sex 
Classification 



190 



191 



192 




194 



FIG. 19 



SUBSTITUTE SHEET (RULE 26) 



INTERNATIONAL SEARCH REPORT 



Internatiorfai/lppHcation No 

PCT/US 00/34780 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 7 A61B5/00 G01N21/35 



According to International Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 7 A61B G01N 



Documentation searched other than minimum documentation to the extent that such documents are included In the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 

WPI Data, PAJ, EPO-Internal 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 0 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



WO 89 08428 A (FUTREX INC) 
21 September 1989 (1989-09-21) 



claims 1,2,7; figure 1 

US 5 013 906 A (MIYAKAWA TATSU0 ET AL) 

7 May 1991 (1991-05-07) 

cited In the application 

column 2, line 62 -column 3, line 29 

-/- 



1-5,14, 
21,22, 
26,27, 
. 30,32,33 



1,14,21, 

26,30, 

32,33 



II 



Further documents are listed in the continuation of box C. 



ID 



Patent family members are listed in annex. 



° Special categories of cited documents : 

"A" document defining the general state of the art which is not 
considered to be of particular relevance 

•E' earlier document but published on or after the International 
filing date 

■L" document which may throw doubts on priority claim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

•O" document referring to an oral disclosure, use, exhibition or 
other means 

■p" document published prior to the International filing date but 
later than the priority date claimed 



T later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
Invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

■Y" document of particular relevance; the claimed Invention 
cannot be considered to Involve an Inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

document member of the same patent family 



Date of the actual completion of the IntemalionaJ search 



27 April 2001 



Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan 2 
NL-2280HVRijswijk 
TeL (431-70) 340-2040, Tx. 31 651 eponl, 
Fax: (+31-70) 340-3016 



Date of mailing of the International search report 



09/05/2001 



Authorized officer 



Krametz, E 



Foim PCT/ISA/210 (second sheet) {July 1992) 



page 1 of 2 



INTERNATIONAL SEARCH REPORT Internatio Plication No 

PCT/US 00/34780 


C(Contlnuatlon) DOCUMENTS CONSIDERED TO BE RELEVANT 


Category 0 


Citation of document, with Indication .where appropriate, of the relevant passages 


Relevant to claim No. I 


A 


SHUSTER S ET AL: "THE INFLUENCE OF AGE 
AND SEX ON SKIN THICKNESS, SKIN COLLAGEN 
AND DENSITY" 

BRITISH JOURNAL OF DERMATOLOGY . , 

vol. 93, 1975, pages 639-643, XP000998738 

ISSN: 0007-0963 

cited in the application 

page 642, paragraph 2 - paragraph 4 


1,4,5, 
14,21, 
26,30, 
32,33 



Form FCT/ISA/210 (continuation of second sheet) (July 1992) 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 

Information on patent family members 



Patent document 
cited in search report 



Publication 
date 



mternatio^ai^pplicatlon No 

PCT/US 00/34780 



WO 8908428 



21-09-1989 



US 5013906 



07-05-1991 



Patent family 
member(s) 



Publication 
data 



US 


4850365 A 


25- 


-07- 


1989 


us 


4928014 A 


22- 


-05- 


1990 


AT 
A! 




1 R- 






Al 1 
AU 




UD 




.IQftQ 


DE 


68927321 D 


14- 


•11- 


•1996 


DE 


68927321 T 


13- 


•02- 


■1997 


EP 


0408637 A 


23- 


-01- 


•1991 


JP 


2648377 B 


27- 


-08- 


•1997 


JP 


4500762 T 


13- 


-02- 


■1992 


us 


4990772 A 


05- 


-02- 


■1991 


JP 


2076528 A 


15- 


-03- 


-1990 


CA 


1306782 A 


25- 


-08- 


■1992 



Rxm PCT/ISA/210 (patent family arm ox) (July 19S2) 



WO 01/52726 PCT/USOO/34780 

In one embodiment a group of LEDs are employed to transmit energy at pre- 
selected wavelengths to the skin and radially surround a single detection 
element at specific distances. The LEDs are alternately energized and the 
detected energy of each LED that is reflected or transmitted through the skin 
5 is used to form one spectrum. While numerous distances are possible, the 
preferred implementation has a minimum of 1 mm and maximum of 3 mm 
edge-to-edge distance between the LEDs and the point of detection. The set 
of wavelengths include 1070, 1220, 1250, 1330, 1380, 1430 nanometers 
and/or 1550, 1560, 1670 and 1800 nanometers. Coupling of the illumination 
10 and detection elements is accomplished through starring optics and a lens 
system. One skilled in the art can appreciate that other coupling methods are 
also applicable including fiber optics given the criterion for the distance 
between the point of illumination and detection. 

1 5 In the preferred embodiment of the invention, the instrument employs a quartz 
halogen lamp, a monochrometer, and InGaAs detectors. The detected 
intensity from the sample is converted to a voltage through analog electronics 
and digitized through a 16-bit A/D converter 17. The spectrum is passed to 
the sex classification procedure 100 for processing. First, the absorbance is 

20 calculated (101) on the basis of the detected light through -\og(R/R 0 ) where R 
is the reflected light and R 0 is the light incident on the sample determined by 
scanning a reference standard. Subsequent processing steps, described 
below, result in either a sex determination or a message indicating an invalid 
scan. 

25 

Alternately, the measurement can be accomplished with existing NIR 
spectrometers that are commercially available, including a Perstorp Analytical 
NIRS 5000 spectrometer or a Nicolet Magna-IR 760 spectrometer. In 
addition, the measurement can be made by collecting reflected light off the 
30 surface of the skin or light transmitted through a portion of the skin, such as 
the finger or the ear lobe. Further, the use of reflectance or transmittance can 
replace the preferred absorbance measurement. 
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